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Abstract 

The minimal length of "one-dimensional" Josephson junctions, in which the specific bound 
states of the magnetic flux retain their stability is discussed numerically. Thereby, we consider as 
"long" every Josephson junction, in which there exists at least one nontrivial stable distribution 
of the magnetic flux for fixed values of all the physical and the geometrical parameters. 

Our results can be applied for optimization of the sizes of devices containing JJs for different 
operating conditions. 

1 Formulation of the Problem 

It is known that the stationary distributions (bound states) of the magnetic flux ip(x) in "long" (one- 
dimcnsional) Josephson junctions (JJs) can be modeled as solutions of the nonlinear boundary value 
problem (BVP): 

- (p xx + jD{x)smip + = 0, ip x (±R) = h B , (1) 

where 7 is the bias current, and hs is the external (boundary) magnetic field. We note that the kind 
of boundary conditions, either in presence or in absence of the current 7 in the right-hand side of eqn. 

t, is determined by the geometry of the junction. Here we consider simple JJs with overlap geometry 
, in which the current 7 can be approximately considered as a constant. The generalization of our 
results in the case of another geometry, for example in-line geometry, does not require big efforts. 

We suppose that the given function < jc(x) < 1 describes the variations of the Josephson 
current amplitude, caused by the possible local inhomogeneities of the dielectric layer thickness p[. 
When Jd(x) = 1 the junction is homogeneous (see the complete analytical results in recent paper B). 

The single inhomogeneity is characterized by its respective position xo and width fi, so the amplitude 
jo = 3d(x, xq, fj,). The frequently used expression jr> (x) = 1 — [j,5(x— xq) (here S(x) is Dirac's function) 
is useful in theoretical examination (see the basic papers Q and ||). In this paper we use an isosceles 
trapezium with base jj, (Fig.l) as a more suitable physical model of inhomogeneity M. fj\. 
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Figure 1 



The stability of some concrete solution (p(x) of the BVP is determined by the sign of the minimal 
eigenvalue X m in of the corresponding Sturm-Liouville problem (SLP): 



ipxx + q(x) ip = X ip , ij) x (±R)=0, 



(2) 



where q(x) = Jd{x) cos y(x) is the potential of SLP, originated from the mentioned solution <p(x). 

Each of the different eigenvalues \ m in = Ao < Ai < A2 < . . . A n < . . . , conforms to an unique 
eigenfunction ip„(x) , n = 0, 1, 2, . . . , which satisfies some additional (norm) condition, for example: 



ip 2 (x) dx = 1 . 



(3) 



If the minimal eigenvalue X rn in > , then the respective solution ip(x) of BVP (0) is stable with 
respect to small time-space perturbations. If the minimal eigenvalue X m in < 0, then this solution is 
unstable. The value X m in — corresponds to a bifurcation point, in which the stable solutions of eqn. 
(|l|) go to unstable ones and vice versa (for details see ||). 

Apart from the space coordinate x, the virtual solutions of the nonlinear BVP (|l|) depend also 
on the physical parameters hs, 7 and the "technological" ones /1, A = 2R (length of the junction), 
and Xq, i.e., tp = ip[x,p), where we simply substitute p = {hg, 7, /i, xq, A}. The varying of each of 
those parameters causes a variation of the distribution (p(x,p) and, therefore, subsequent variations 
of the potential q(x,p), the eigenvalues A(p), and the respective eigenfunctions ip(x,p). Thus, we can 
conclude that every concrete solution of BVP (|J) has an area, where this solution remains stable with 
regard to the variations of the parameters p. In the parametric space the equation: 



Xmin ip) = 



(4) 



determines a hypersurface, whose points appear to be bifurcation points corresponding to the solution 
under consideration. The intersections of the bifurcation hypersurface (||), when there are chosen pairs 
of the parameters p for fixed values of others, are called bifurcation curves and the respective values of 
the parameters - bifurcation (critical) parameters. From a physical viewpoint most interesting seem to 
be the bifurcation curves of the kind "current - magnetic flux" A m j n (/ifl,7) = 0, when the geometrical 
parameters /i, R, and xq are fixed. These curves could be obtained experimentally [H 0, [Tof. 
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Case: A = 6.3, (x = 1, h B = 0, y = 




Length of the Junction 

Figure 2 

From a technological point of view, however, it is worth investigating the bifurcation curves as 
functions with respect to at least one of the geometrical parameters, for example: 

X min (h B ,A) = 0,A min (7, A) = , X min ([i, A) =0, ... 

Such kind of problems can be connected with the optimization of sizes of devices, containing JJs. 

Most of our results below are related to the typical solution of BVP ([!]) , physically existing both 
in homogeneous and inhomogencous JJs: so called single fluxon. Let us note that in the infinite 
homogeneous or inhomogeneous JJ with one <5-shaped microinhomogeneity at the point xq — 0, the 
single fluxon/antifluxon is represented by the exact solution [||, jjj: 

(f(x) — 4 arctanexp(±x). (5) 

More complicated ( mult i- fluxon) solutions there exist in inhomogeneous JJs for big enough values 
of He- 

The numerically obtained basic functions: the magnetic flux <p(x), the magnetic field ip x (x), and 
the first eigenfunction ipo(x) of such solution for A = 6.3, /i = 1, hs = 0, and 7 = are represented 
at Fig. 2. 

Fig. 3 represents the typical relationships A m j n (A), corresponding to the single fluxons in inho- 
mogeneous junctions with one resistive inhomogencity placed at the point xq = when /i = 0.5 and 
[i=l, respectively. The points Bq and Bi appear to be bifurcation points for these solutions, i.e., 
the zeroes of the corresponding relationships A m i„(A). According to the above mentioned reasonings, 
these points determine the minimal length of the J J, for which the single fluxon is still stable for given 
values of the other parameters. 
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Usually, from physical point of view, the concrete JJ is "long" or "one-dimensional" if the Josephson 
penetration depth Xj satisfies the condition Xj >> 1. However, such definition is not perfect enough 
mathematically. On the other hand, it seems that the calculation of the minimal geometrical sizes and, 
especially, the minimal length A of homogeneous or inhomogeneous JJs, corresponding to a concrete 
nontrivial distribution of the magnetic flux, is an important practical problem. 

The trivial method, obviously, is through the point by point construction of the curves of type 
A m i„(A) using BVPs (|l|)-(^). However, it is quite natural to put the question how to calculate directly 
the bifurcation curves. A general approach of study for the formulated problem is proposed in the 
papers @ p. 

We consider the system of the equations ([!]) , (|^) , and (^) as a closed nonlinear eigenvalue problem 
with respect to the functions ip(x), ip(x) and one of the parameters p, for example, hs or 7 in ]l2| , 
or R in this case, while the other parameters, and the eigenvalue A as well, are given. Fixing A > 
to be small enough (for example A = 0.001 -j- 0.01), we guarantee a stable solution to be obtained. 
We note, that the derivative dX/dpo — > 00 when po approaches its critical (bifurcation) value, so, 
the solutions of the above system with a priori prescribed accuracy belong to the small vicinity of the 
sought bifurcation curve. 

The generalized Continuous Analog of Newton's Method jL3) for solving the mentioned system is 
used (see details in our recent paper 0). 

We have to note that for some initial approximations the method can lead to an unstable solution, 
as well. In this case, the given A corresponds to one of the higher eigenvalues, not to the minimal one. 
The simplest method to avoid this difficulty uses the oscillation theorem (m) for eigenfunctions: // the 
obtained function ip( x ) is the minimal eigenfunction, then it does not have a zero on the interval A. 



2 Results and Discussion 

Fig. 4 shows the convergence of the iteration process. The upper curve, marked by "V" , presents the 
initial distribution of the magnetic field ip x {x) alongside the junction. This distribution is chosen to 
be a solution of BVP (|l|), for A = 10, /i = 1, hs = 0, and 7 = 0. At the same figure the lower curve, 
marked by "o" , presents the calculated distribution of the magnetic field corresponding to the minimal 
eigenvalue X m i n = 0.01. Hence, this is the sought bifurcation distribution and, in the framework of 
this model, the value of the spectral parameter A m i n 4.24 is the minimal length of the junction, 
providing a stable single fluxon. In this sense, junctions, whose length lies below the critical value 
A mi „ of A, have to be considered as "short" for distributions <p(x) of the kind single fluxon. However, 
for such a short length A there are a unique stable and unstable distributions of the magnetic flux in 
the junction - Meissner's solutions: the trivial solutions <p(x) = 0, 2ir, . . . (stable) and tp(x) — n, 3tt, . . . 
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(unstable) for Jib = and 7 = 0. 



Case: h B = 0, y = 0, (i = 1 




Length of the Junction 



Figure 4 

The influence of the parameters Iib, 7, and xo on the minimal length A m j„ of the JJs will be 
considered below. 

If the current 7^0 and the values of all other parameters are fixed, then the stability of the fluxon 
requires the junction to be longer (see Fig. 5). This fact corresponds to the classical conclusion pj that 
the existence of quantity 7 violates the stability of the bound states in JJs. The "force" 7 displaces 
the maximum of ip x (x) from Xq = to the left or right depending on its sign; the full magnetic flux 
Aip = <p(R) — (p(—R) increases, so that the minimal junction's length A m j„, necessary to accommodate 
such a flux, increases as well. 



Case: h B = 0, ja = 1 




Length of the Junction 



Figure 5 



5 



Case: /z R = 




0.06 0.12 0.18 0.24 

Current y 

Figure 6 

The numerically obtained influence of the current 7 upon the bifurcation length A min is plotted on 
Fig. 6 for two values of the size fj, of the inhomogeneity: (i = 0.5 and ji = 1. The curves increase very 
abruptly and when 7 approaches some critical value, the derivative dA min /d"f — ^ 00. This behaviour 
means that beyond the critical value of 7, JJs which correspond to the single fluxon do not exist. 



Case: y = 




0.5 1 1.5 2 2.5 

Boundary Magnetic Field h B 

Figure 7 

On the contrary, if the boundary magnetic field hs ^ and varies within some fixed domain, then 
the stability of the fluxon is improved (see Fig. 7). Well noticeable is the stabilizing influence of the 
boundary magnetic field hs upon the critical length of the junction, i.e., the increase of the magnitude 
of hs decreases the minimal length of the junction A m i„. Thus, the calculated minimal length of the 
J J corresponding to the single fluxon = 0, 7 = 0) is A m i n ~ 4.24, while the same quantity when 
hs = 1 and 7 = is A. m i n ~ 2.7. 
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Fig. 8 represents the distributions of magnetic field f x (x) alongside the junction in absence of 
current (7 = 0) depending on its boundary values hg = and hg = 1, respectively. It is clear that 
the two solutions could be well approximated by means of two-degree polynomials. Furthermore, they 
seem to be geometrically similar. Namely, the curve <-p x {x) corresponding to a magnetic field hg > 
may be considered as obtained from the other curve by a homothety alongside x axis with respect to 
the pole x = and a translation alongside the axis ip x . 



Case: Homogeneous Junction, y = 
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Figure 9 



At the next Fig. 9, the dependence of the minimal length A min on the boundary magnetic field hg 
for 7 = in the case of a homogeneous J J (the amplitude ji>(x) = 1) is presented. Even though this 
dependence is qualitatively similar with the analogous ones at Fig. 7, it differs from the big slope of 
Amin for small hs, so that A m j n (/ie) — > 00 when hg — > 0. 

Such behaviour of the curve A m ,- n (/iB) can be qualitatively explained by means of the formula H: 



^0 — ^min — 



1 
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(G) 



which is exact for the single fluxon (g) in an infinite JJ with one resistive (5-shaped inhomogeneity at 
the center x = 0. Hence, for /i = we obtain X m i n — 0, so that for hg — and 7 = the single 
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fluxon in homogeneous JJ has to be considered as a quasi-stable (bifurcation) state of the magnetic 
flux ip(x) ||. On the contrary, in inhomogeneous J J with resistive inhomogeneities the single fluxon is 
still stable (\ m in > 0), which physically is due to the existence of the inhomogeneity. 

In the case of an inhomogeneous J J with finite length A = 6.3 the calculated curve X m i„(fi) (marked 
by "□") for hs — and 7 = is shown at Fig. 10. We note that for small enough values of the size /x 
the dependence is analogous to (|J). 

Let us introduce the full energy F connected with some solution ip(x) of BVP (fil) via the formula: 



;vl + 3d (x) (1 - cos ip) + ~/tp 



dx + h B [<p{R)-cp{-R)}. 



It is efficient to compare the energy curves F(fj,) (marked by "V") for fixed A = 6.3 and F m i n (n) 
(marked by "o") for a JJ with minimal length (see Fig. 10). The magnitudes of the energy are 
regarded to the energy (F = 8) of the single fluxon (g) in an infinite homogeneous JJ. It is well 
noticeable that for /i > the inequality F(fi) > F min (^i) is satisfied, e.g., the single fluxon in the JJ 
with minimal length is more stable. The cusp for fi = corresponds to a bufurcation in accordance 
with formula (n). 



Case: h B = 0, y = 
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Width of the Inhomogeneity ]x 



Figure 10 



Case: h B = 0, y = 




0.4 0.8 1.2 1.6 2 

Width of the Inhomogeneity \x 

Figure 11 
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At Fig. 11 the numerically obtained relationship A TO j„(//) is drawn, depending on the length of the 
trapezium's base y, (width of the inhomogeneity) (see Fig.l). The presented three curves corresponding 
to different locations of the centre of inhomogeneity Xo are qualitatively similar. The most important 
their feature is the existence of minimums, which are placed in the interval fi £ (0.9, 1.3). Therefore, 
there exist such values of the quantity fj,, which provide a minimal bifurcation length of a JJ from 
among all the minimal lengths for given other parameters. Let us call them "optimal" lengths. For 
example, if xq — ("symmetric" inhomogeneous JJ), the optimal length A opt = A min ~ 4.22 of the 
JJ corresponds to [i ~ 1.1, whereas A m i„ ~ 4.48 if the width of the inhomogeneity is /i = 0.5, and 
A mm ~ 4.3 for /z = 1.5. 



Case: h B = 0, jj, = 1 




-0.8 -0.4 0.4 0.8 

Coordinate of Inhomogeneity x Q 



Figure 12 

The bifurcations of the bound states in the circular JJs with respect to the location xq of the 
(5-shaped inhomogeneity are investigated numerically in jl^]. Fig. 12 shows the curves A m i n (xo) de- 
scribing the behaviour of the minimal length of the linear JJs. When the current 7 = and the width 
of inhomogeneity ji = 1 is fixed, the curve is symmetric and, apparently, the minimal value of A m i n is 
reached at the location xo — 0. When the current becomes nontrivial (7 = 0.1), the minimum shifts 
into some other point. For negative values of 7 the dependence is the same but mirror reflected with 
respect to the line xq = 0. 



Case: \i = 1, h n = 3.2, y = 
3.8 -i 
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Figure 13 
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Finally, the last Fig. 13 provides as an illustration a three-soliton bifurcation solution. Such kind 
of solutions can be obtained when the boundary magnetic field hs is big enough. Here, the concrete 
solution is derived for hs — 3.2. At the same graph the initial distribution of the magnetic field ip x (%) 
as solution of BVP ([!]), for A = 8, 7 = 0, and fx = 1 is included. 
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